Scanamorphos: a map-making software for Herschel and similar scanning bolometer arrays

Scanamorphos is one of the public softwares available to post-process scan observations performed with the Herschel photometer arrays. This post-processing mainly consists in subtracting the total low-frequency noise (both its thermal and non-thermal components), masking high-frequency artefacts such as cosmic ray hits, and projecting the data onto a map. Although it was developed for Herschel, it is also applicable with minimal adjustment to scan observations made with some other imaging arrays subjected to low-frequency noise, provided they entail sufficient redundancy; it was successfully applied to P-Artemis, an instrument operating on the APEX telescope. Contrary to matrix-inversion softwares and high-pass filters, Scanamorphos does not assume any particular noise model, and does not apply any Fourier-space filtering to the data, but is an empirical tool using purely the redundancy built in the observations -- taking advantage of the fact that each portion of the sky is sampled at multiple times by multiple bolometers. It is an interactive software in the sense that the user is allowed to optionally visualize and control results at each intermediate step, but the processing is fully automated. This paper describes the principles and algorithm of Scanamorphos and presents several examples of application.Comment: This is the final version as accepted by PASP (on July 27, 2013). A copy with much better-quality figures is available on http://www2.iap.fr/users/roussel/herschel

Recent Development of MESSINE, a 3D Eddy Current Model

In the 1996 QNDE Conference, we presented a parametric forward model [1], which has been recently named MESSINE (Model for Electromagnetic Simplified Simulation In Nondestructive Evaluation), to predict eddy current signal. The proposed model first discretizes the eddy current distribution into current loops. A parametric description of the shape of these loops is given according to the observation of the results obtained with a three-dimensional finite element code which provides a realistic distribution of the induced currents. The loops’ inductances and resistances are then calculated. By considering the system constituted of the coil and the current loops as a « multi-transformer», their current intensity is determined. The impedance change, which is the component of the eddy current signal, can then be deduced. The model was validated in the case of axisymmetric configurations. Comparisons with both analytical (Dodd and Deeds [2]) and numerical models showed very good agreements. Then the proposed model was applied to three-dimensional configurations. Impedance changes of a coil along rectangular through-wall slot were calculated. Comparisons with experimental results show a fairly good agreement for the impedance change phases, but a poorer one for the impedance change amplitudes. Investigations were made to improve the parametric description of the current loop deformation. One of the solutions to improve the parametric description is presented here

Evidence of triggered star formation in G327.3-0.6. Dust-continuum mapping of an infrared dark cloud with P-ArT\'eMiS

Aims. Expanding HII regions and propagating shocks are common in the environment of young high-mass star-forming complexes. They can compress a pre-existing molecular cloud and trigger the formation of dense cores. We investigate whether these phenomena can explain the formation of high-mass protostars within an infrared dark cloud located at the position of G327.3-0.6 in the Galactic plane, in between two large infrared bubbles and two HII regions. Methods: The region of G327.3-0.6 was imaged at 450 ? m with the CEA P-ArT\'eMiS bolometer array on the Atacama Pathfinder EXperiment telescope in Chile. APEX/LABOCA and APEX-2A, and Spitzer/IRAC and MIPS archives data were used in this study. Results: Ten massive cores were detected in the P-ArT\'eMiS image, embedded within the infrared dark cloud seen in absorption at both 8 and 24 ?m. Their luminosities and masses indicate that they form high-mass stars. The kinematical study of the region suggests that the infrared bubbles expand toward the infrared dark cloud. Conclusions: Under the influence of expanding bubbles, star formation occurs in the infrared dark areas at the border of HII regions and infrared bubbles.Comment: 4 page

Probing the structure of a massive filament: ArTéMiS 350 and 450 μm mapping of the integral-shaped filament in Orion A

Context. The Orion molecular cloud is the closest region of high-mass star formation. It is an ideal target for investigating the detailed structure of massive star-forming filaments at high resolution and the relevance of the filament paradigm for the earliest stages of intermediate- to high-mass star formation. Aims. Within the Orion A molecular cloud, the integral-shaped filament (ISF) is a prominent, degree-long structure of dense gas and dust with clear signs of recent and ongoing high-mass star formation. Our aim is to characterise the structure of this massive filament at moderately high angular resolution (8′′ or ~0.016 pc) in order to measure the intrinsic width of the main filament, down to scales well below 0.1 pc, which has been identified as the characteristic width of filaments. Methods. We used the ArTéMiS bolometer camera at APEX to map a ~0.6 × 0.2 deg2 region covering OMC-1, OMC-2, and OMC-3 at 350 and 450 μm. We combined these data with Herschel-SPIRE maps to recover extended emission. The combined Herschel-ArTéMiS maps provide details on the distribution of dense cold material, with a high spatial dynamic range, from our 8′′ resolution up to the transverse angular size of the map, ~10-15′. By combining Herschel and ArTéMiS data at 160, 250, 350, and 450 μm, we constructed high-resolution temperature and H2 column density maps. We extracted radial intensity profiles from the column density map in several representative portions of the ISF, which we fitted with Gaussian and Plummer models to derive their intrinsic widths. We also compared the distribution of material traced by ArTéMiS with that seen in the higher-density tracer N2H+(1-0) that was recently observed with the ALMA interferometer. Results. All the radial profiles that we extracted show a clear deviation from a Gaussian, with evidence for an inner plateau that had not previously been seen clearly using Herschel-only data. We measure intrinsic half-power widths in the range 0.06-0.11 pc. This is significantly larger than the Gaussian widths measured for fibres seen in N2H+, which probably only traces the dense innermost regions of the large-scale filament. These half-power widths are within a factor of two of the value of ~0.1 pc found for a large sample of nearby filaments in various low-mass star-forming regions, which tends to indicate that the physical conditions governing the fragmentation of pre-stellar cores within transcritical or supercritical filaments are the same over a large range of masses per unit length. © F. Schuller et al. 2021

Characterizing filaments in regions of high-mass star formation: High-resolution submilimeter imaging of the massive star-forming complex NGC 6334 with ArTeMiS

Context. Herschel observations of nearby molecular clouds suggest that interstellar filaments and prestellar cores represent two fundamental steps in the star formation process. The observations support a picture of low-mass star formation according to which filaments of ~0.1 pc width form first in the cold interstellar medium, probably as a result of large-scale compression of interstellar matter by supersonic turbulent flows, and then prestellar cores arise from gravitational fragmentation of the densest filaments. Whether this scenario also applies to regions of high-mass star formation is an open question, in part because the resolution of Herschel is insufficient to resolve the inner width of filaments in the nearest regions of massive star formation. Aims. In an effort to characterize the inner width of filaments in high-mass star-forming regions, we imaged the central part of the NGC 6334 complex at a resolution higher by a factor of >3 than Herschel at 350 μm. Methods. We used the large-format bolometer camera ArTéMiS on the APEX telescope and combined the high-resolution ArTéMiS data at 350 μm with Herschel/HOBYS data at 70–500 μm to ensure good sensitivity to a broad range of spatial scales. This allowed us to study the structure of the main narrow filament of the complex with a resolution of 8″ or <0.07 pc at d ~ 1.7 kpc. Results. Our study confirms that this filament is a very dense, massive linear structure with a line mass ranging from ~500 M⊙/pc to ~2000 M⊙/pc over nearly 10 pc. It also demonstrates for the first time that its inner width remains as narrow as W ~ 0.15 ± 0.05 pc all along the filament length, within a factor of <2 of the characteristic 0.1 pc value found with Herschel for lower-mass filaments in the Gould Belt. Conclusions. While it is not completely clear whether the NGC 6334 filament will form massive stars in the future, it is two to three orders of magnitude denser than the majority of filaments observed in Gould Belt clouds, and has a very similar inner width. This points to a common physical mechanism for setting the filament width and suggests that some important structural properties of nearby clouds also hold in high-mass star-forming regions

Measurements of the Suprathermal Electron Density Profile by Electron Cyclotron Emission at the Upper Cut-Off Layer in Tore Supra

The possibility of measuring the radial profile of fast-electron distributions from ECE spectra is discussed. The most usual ECE diagnostic set-up, i.e. the second-harmonic X-mode, observed by low-held side antennae in the equatorial plane, can be used by exploiting the reflection of spontaneously emitted radiation by the upper cut-off layer and the length of the sightline between the reflecting surface and the receiving antenna. The emitted intensity depends on this length, which introduces an element of spatial resolution in the suprathermal ECE spectra. Two methods are proposed. First, a dynamic method is described in which the length of the sightline is varied by density perturbations. Second, a spectral method is discussed that uses the length variation of the sightline with the observation frequency. A proof of the principle behind both methods is given for a lower hybrid current driven discharge in Tore Supra in which a pellet was injected

Recent Development of MESSINE, a 3D Eddy Current Model

In the 1996 QNDE Conference, we presented a parametric forward model [1], which has been recently named MESSINE (Model for Electromagnetic Simplified Simulation In Nondestructive Evaluation), to predict eddy current signal. The proposed model first discretizes the eddy current distribution into current loops. A parametric description of the shape of these loops is given according to the observation of the results obtained with a three-dimensional finite element code which provides a realistic distribution of the induced currents. The loops’ inductances and resistances are then calculated. By considering the system constituted of the coil and the current loops as a « multi-transformer», their current intensity is determined. The impedance change, which is the component of the eddy current signal, can then be deduced. The model was validated in the case of axisymmetric configurations. Comparisons with both analytical (Dodd and Deeds [2]) and numerical models showed very good agreements. Then the proposed model was applied to three-dimensional configurations. Impedance changes of a coil along rectangular through-wall slot were calculated. Comparisons with experimental results show a fairly good agreement for the impedance change phases, but a poorer one for the impedance change amplitudes. Investigations were made to improve the parametric description of the current loop deformation. One of the solutions to improve the parametric description is presented here.</p